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^p ' The resolution of the problem of cosmological singularity in the frame- 

^ , work of gauge theories of gravitation is discussed. Generalized cosmo- 

^^1 logical Friedmann equations for homogeneous isotropic models filled by 

interacting scalar fields and usual gravitating matter are deduced. It 
l^£\, is shown that generic feature of cosmological models of fiat, open and 

closed type is their regular bouncing character. 



$H ' 1.1. Introduction 



As it is well known the problem of cosmological singularity (PCS) is one 
of the most principal problems of general relativity theory (GR) and rela- 
tivistic cosmology. There were many attempts to resolve PCS in the frame 
of GR as well as various generalizations of Einstein's theory of gravitation. 
A number of regular cosmological solutions was obtained in the frame of 
metric theories of gravitation and also other theories, in the frame of which 
gravitation is descr ibed by using more general geometry than the Ricman- 
nian one (se d I I I and Refs given herein). In connection with this note 
that the resolution of PCS means not only obtaining regular cosmological 
solutions, but also excluding singular solutions of cosmological equations, 
as result generic feature of cosmological solutions has to be their regular 
character. Moreover, gravitation theory and cosmological equations have 
to satisfy the correspondence principle with Newton's theory of gravitation 
and GR in the case of usual gravitating systems with sufficiently small en- 
ergy densities and weak gravitational fields excluding nonphysical solutions. 
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The greatest part of existent attempts to resolve the PCS does not satisfy 
indicated conditions. Let's take for instance some examples. Although GR 
permits to build regular closed inflationary model s, th e greatest part of 
cosmological inflationary solutions of GR is singulail^'^. Metric theories of 
gravitation based on gravitational Lagrangians including terms quadratic 
in the curvature tensor lead to cosmological equations with high derivatives, 
and although these theories permit to obtain regular cosmological solutions 

7 

with Friedmann and de Sitter asymptotics , however, they possess non- 
physical solutions also. At last time single regular cosmological solutions 
were found in the frame of superstring theory (brane cosmology) by us- 
ing speciflc suppositions (see for example*"^^). To build cyclic model of the 
Universe specific negative scalar field potential was introduced in Ref j I I 

As it was shown in a number of our papers^^"^**' ^^^, the gauge theories 
of gravitation (GTG) - the Poincare GTG, metric-affine GTG - possess 
important regularizing properties and permit to resolve the PCS. Satisfy- 
ing the correspondence principle with GR in the case of usual gravitat- 
ing systems with rather small energy densities and pressures, the GTG 
lead to essentially different consequences in the case of gravitating systems 
at extreme conditions with extremely high energy densities and pressures, 
namely the GTG lead to conclusion on possible existence of limiting energy 
density for gravitating systems, nearby of which gravitation has the char- 
acter of repulsion, but not attractiorP^. The investigation of homogeneous 
isotropic cosmological models in GTG shows, that by certain restrictions on 
equation of state of gravitating matter at extreme conditions and indefinite 
parameter of cosmological equations of GTG generic feature of cosmologi- 
cal models including infiationary models is their regular bouncing character. 
Such conclusions were obtained, in particular, in our recent papers *' for 
inflationary cosmological models filled by noninteracting scalar fields and 
gravitating matter with linear equation of state. 

Present paper is devoted to further study of regular cosmology in 
the frame of GTG. In Section 2 cosmological equations for homogeneous 
isotropic models filled by interacting scalar fields and gravitating matter 
with equation of state in general form are deduced. In Section 3 general 
mathematical properties of deduced equations and their solutions are an- 
alyzed. As illustration of discussed theory in Section 4 particular solution 
for regular inflationary cosmological model is given. 
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1.2. Generalized cosmological Friedmann equations in GTG 

Homogeneous isotropic models in GTG are described by tlie following gen- 
eralized cosmological Friedmann equations (GCFE) 



k r d , 
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where R{t) is the scale factor of Robertson- Walker metrics, k ~ +1,0,-1 
for closed, flat, open models respectively, /3 is indefinite parameter with 
inverse dimension of energy density, Mp is Planckian mass, a dot denotes 
differentiation with respect to time''. (The system of uni ts w ith S = c = 1 is 
used) . At first the GCFE were deduced in Poincare GTGM^, a nd late r it was 
shown that Eqs.(l)-(2) take place also in mctric-afhne GTGpSESI. in the 
frame of considered theory the conservation law in usual form takes place 

p + 3H{p + p)^0, (3) 

where H = ^ is the Hubble parameter. Besides cosmological equations 
(l)-(2) gravitational equations of GTG lead to the following relation for 
torsion function S and nonmetricity function Q 

^-^Q = -J|ln|l-/3(p-3p)|. (4) 

In Poincare GTG Q = and Eq. (4) determines the torsion function. 
In metric-affine GTG there are three kinds of models^-* in the Riemann- 
Cartan space-time ((5 = 0), in the Weyl space-time (5 = 0), in the Weyl- 
Cartan space-time {S ^ 0, Q ^ 0, the function S is proportional to the func- 
tion Q). The value of \(3\~^ determines the scale of extremely high energy 
densities. The GCFE (l)-(2) coincide practically with Friedmann cosmolog- 
ical equations of GR if the energy density is small \(3{p — 3p) | <C 1 (p ^ ^p). 
The difference between GR and GTG can be essential at extremely high 
energy densities \/3{p — 3p)\ > 1. Ultrarelativistic matter {p = ^p) and 
gravitating vacuum {p = — p) with constant energy density are exceptional 
systems because Eqs. (l)-(2) are identical to Friedmann cosmological equa- 
tions of GR in these cases independently on values of energy density and 
S = Q ~ 0. The behaviour of solutions of Eqs. (l)-(2) depends essentially 



^Parameter f3 is defined as /3 = — ^ (IOtt) f Mp , where / is linear combination of 
coefficients at terms of gravitational Lagrangian quadratic in the curvature tensor. 
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on equation of state of gravitating matter at extreme conditions and on sign 
of parameter (3. In the case of models filled by gravitating matter without 
scalar fields the GCFE (l)-(2) lead to regular in metrics cosmological so- 
lutions by the following restrictions: 1) /3 > and at extreme conditions 
p < |p, 2) /? < and p > ^p^. The investigation of models including scalar 
fields on the base of Eqs (l)-(2) shows that the choice /3 < permits to 
exclude the divergence of time derivative of scalar fields^. In connection 
with this we put below that parameter /3 is negative and in the frame of 
our classical description \P\~^ < 1 • M^. 

By using GCFE (l)-(2) we will study homogeneous isotropic models 
filled by interacting scalar field </> minimally coupled with gravitation and 
gravitating matter with equation of state in general form pm = Pm{Pm)^- 
Then the energy density p and pressure p take the form 



1 



i2 



p = ^r + v + pn 



1 -9 

P=7;<P -V + Pr^ 



(5) 



where scalar field potential V — V{4),pm) includes the interaction between 
scalar field and gravitating matter. In the most important particular case 
of radiation (ultrarelativistic matter) the expressions of 1^(0, pm) can be 
obtained by taking into account temperature corrections for given scalar 
field potentials ' and the following relation for energy density pm ^ T^ {T 
is absolute temperature) . By using the scalar field equation in homogeneous 
isotropic space 



-*-§^ 



(6) 



we obtain from Eqs. (3), (5), (6) the conservation law for gravitating matter 



dV 
dp,n 



+ 3H {pm +Pm) = 0. 



(7) 



By virtue of Eqs.(5)-(7) the GCFE (l)-(2) can be transformed to the fol- 
lowing form 







{3pm ~W)-1 



3/3^0 1 



i?2 3M2 
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Z (8) 



The generalization for the case with several scalar fields can be made directly. 
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where Z — 1 — (3 i AV — (fP' + p„i — 3pm ) • Relation (4) takes the form 



S--Q 



3/3 
2Z 



H 



i2 



I P 



2 1 



ni I I-'ni 



dV 
9p™ 



rfpr, 



(3p™ - 4y) - 1 



dV 



(10) 

Unlike GR the cosmological equation (8) leads to essential restrictions on 
admissible values of scalar field and gravitating matter. Imposing /3 < 0, 
we obtain from Eq. (8) in the case fc = 0, +1 



Z >0, 



<AV 



+ Pm- 3Pn 



(11) 



The condition (11) is valid also for open models discussed below. This means 
that relation (11) is fulfilled for all cosmological models independently on 
their type (fc — 0, +1,— 1). Now let us introduce the 3-dimensional space 
P with axes {(f>,(j), pm)- The domain of admissible values of scalar field 0, 
time derivative (j) and energy density pm in space P determined by (11) is 
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limited by bound L defined as 

Z = or ^ ^ ± {W + \I3\-^ + prr, - 3pmy . (12) 

From Eq. (8) the Hubble parameter on the bound L is equal to 

dV- 

Hl = . „ , . .^ , ^ r- (13) 



1 Pn 



2 1 



/ d 



dV \ du 

dp:: ^ ^" 



3pn 



AV 



The right-hand part of Eq. (10) with the Hubble parameter determined by 
(13) is equal to \H ^ this means that the torsion (nonmetricity) in this case 
will be regular, if the Hubble parameter is regular. 



1.3. Solutions properties of GCFE 

Let us consider the most important general properties of cosmological solu- 
tions of GCFE (8)~(9). At first note, by given initial conditions for variables 
{4>, (j), pm) and value of R there are two different solutions corresponding to 
two values of the Hubble parameter following from Eq. (8): 



dV ■ r- 



H± 



Z-3/3 
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1 Pr, 
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where 



D 



Stt 
3M2 



Pm + l'P^ + V-\(3 (W -(!)''+ p„,- 3p. 



Z 



i?2- 



(15) 
Unlike GR, the values of iJ+ and H- in GTG are sign-variable and, hence, 
both solutions corresponding to if+ and H_ can describe the expansion 
as well as the compression in dependence on their sign. Below we will call 
solutions of GCFE corresponding to if+ and _ff_ as iJ+-solutions and iJ_- 
solutions respectively. Any cosmological solution contains both _ff_- and 
iJ+-solution. In points of bound L we have D = 0, H+ = H- and the Hub- 
ble parameter is determined by (13). The GCFE (8)-(9) are satisfied on the 
bound L, corresponding solutions of GCFE - L-solutions - are their par- 
ticular solutions. The scalar field (p and energy density pm for L-solutions 
can be found by integration of Eqs. (6)^(7), where the Hubble parameter is 
defined by (13), and by using the constrain relation Z = 0. L-solutions can 
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be regular as well singular in dependence on equation of state of gravitat- 
ing matter, scalar field potential V and initial conditions . Trajectories of 
particular solutions situated on the bound L have with 7J±-solutions com- 
mon points, i/_-solutions reach the bound L and Tif 4. -solutions originate 
from them, the surface Z = containing trajectories of particular solutions 
is envelope in space P for cosmological solutions. By using the following 
formula obtained for iJ±-solutions 



7 '^R - 


-f^-- 
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Pm + Pm 


, [iPm 
^dpm V 


w] - 


1] + 2<P^ 




Z-3/3 


^' + 


1 Pm + Pm 
21 + C 


\dpm V 


-41/) 


-1 
/ 





(16) 



it is easy to show that 



lim H± 
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(17) 
From (17) follows that in points of bound L the derivatives i/+ and ij_ 
are equal and their values do not depend on the model type, as result wc 
have the smooth transition from iJ_ -solution to iy_|_-solution on bound i, 
and corresponding cosmological solutions for all types models are regular in 
metrics, Hubble parameter and its time derivative. At the same time accord- 
ing to (17) the value of the time derivative TJ^ for L-solutions is not equal 
to limz~toH± and smooth transition from i7_-solutions to L-solutions and 
from i-solutions to iJ+-sohitions without jump of the derivative H is im- 
possible. Note, that according to Eq. (16) the functions S and Q have the 
following asymptotics for iJ+- and i/_-sohitions at Z ^ 0: 



lim (S Q) — — lim — 



tz- 



(18) 



Unlike flat and open models, for which H^ = H_ only in points of bound 
L and regular inflationary models include H-\.- and iJ_ -solutions reaching 
bound L, in the case of closed models the regular transition from iJ_- 
solution to _ff+-solution is possible without reaching the bound L. It is 
because by certain value of R according to (14)-(15) we have H+ = H- in 
the case Z ^ 0. Such models are regular also in torsion and/or nonmetricity. 
Regular inflationary solution of such type was considered in Ref.l^. 
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In order to study the behaviour of cosmological models at the beginning 
of cosmological expansion, let us analyze extreme points for the scale factor 
R{t): Rq = i?(0), Ho ^ H{0) = 0. (This means that in the case of H+- 
solutions H+o — and in the case of i/_-solutions iJ_o ~ 0). Denoting 
values of quantities at t = by means of index "0" , we obtain from (8)-(9): 



R2 ^0 
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(20) 



where Zq = 1 — /3 I 4Vo — 0o + PmO — ^Pmo ) • A bounce point is described 

by Eq. (19), if the value of Hq is positive. By using Eq.(19) we can rewrite 
the expression of Hq in the form 
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.(21) 



We see from (21) unlike GR the presence of gravitating matter (with 
Pm < Pm) does not prevent from the bounce realization'^ Eq. (19) deter- 
mines in space P extremum surfaces depending on the value of /3 and in the 



'^In GR a bounce is possible only in closed models if the following condition Vq — 4>q 
i(PmO + 3pmo) > takes place. 



February 7, 2008 2:56 



Review Volume: 9in x 6in (Per Chapter) 



regcosml 



Non-Singular Cosmology . . 



case of closed and open models also parametrically on the scale factor Rq. 
In the case of various scalar field potentials applying in inflationary cosmol- 
ogy the value of -ff+o or ^-o is positive on the greatest part of extremum 
surfaces, which can be called "bounce surfaces". '^. By giving concrete form 
of potential V and choosing values of Rq, (po, 4>o and pmo at a bounce, we 
can obtain numerically particular bouncing solutions of GCFE for various 
values of parameter /3. 

The analysis of GCFE shows, that properties of cosmological solutions 
depend essentially on parameter /?, i.e. on the scale of extremely high energy 
densities. From physical point of view interesting results can be obtained, if 
the value of |/3|^^ is much less than the Planckian energy densityl^ i.e. in 
the case of large in module values of parameter /3 (by imposing Mp = 1). In 
order to investigate cosmological solutions at the beginning of cosmological 
expansion in this case, let us consider the GCFE by supposing that 
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Prn ^Pn 



(22) 



Note that the second condition (22) does not exclude that ultrarelativistic 
matter energy density can dominate at a bounce. We obtain: 
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''if 1/3 1 <C Mp the derivative Hq is negative in the neighbourhood of origin of coordi- 
nates in space P that leads to appearance of oscillating solutions of GCFB^sl 
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X 1 
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Equations (23)-(24) do not contain the parameter /3. According to Eq. 
(23) the Hubble parameter in considered approximation is equal to 
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and extreme points of the scale factor are determined by the following 
condition 



Rq 






4Vo - 4>o + PmO - 3pm0 
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(26) 

From Eq. (24) the time derivative of the Hubble parameter at extreme 
points is 
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Obviously Eqs. (26)~(27) correspond to (19)-(20) in considered approxi- 
mation. 

The analysi s giv en in this Section confirms conclusions obtained in our 
previous paperg^l"!. The existence of limiting bound L and bounce surface 
in space P in the case of models including scalar fields and gravitating mat- 
ter ensures regular character of cosmological solutions in metrics, Hubble 
parameter and its time derivative; corresponding restriction on equation 
of state of gravitating matter and potential V at extreme conditions is: 

1 + - — ) ( - — (3p,„ - W) - 1 I > (see (14)). In absence of scalar 
opm J \ dpm J 

fields the bound L is reduced to a point, particular solution of GCFE is 

stationary (_ff=0), and the condition Z = determines a bounce point and 

the value of limiting energy density simultaneously if at extreme conditions 



Pm > 3 A 



„lll 



1.4. Regular inflationary cosmological models 

As illustration of discussed theory regular cosmological models in the sim- 
plest particular case will be considered in this Section. We will consider 
models including noninteractingscalar field with potential V{4>) and ultra- 
relativistic matter {pm = ^Pm] ■ In this case the bound L in space P is 

reduced to two cylindric surfaces (f> — zL (AV + |/?|~^) ^ . Bounce surface is 
reduced also to cylindric surfaces in the case under consideration, when the 
scale of extremely high energy densities is much smaller than the Planckian 
energy density (see (26)). In connection with this we will consider instead of 
space P the plane of variables {(/), (f) and intersections of the bound L and 
bounce surface with this plane. We have in this plane two bound L±-curves 
and in the case of flat models two bounce curves Bi and B2 determined by 
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equation*^ 

Each of two curves Si, 2 contains two parts corresponding to vanishing of 
H+ or iJ_ and denoting by (i?i+, B2+) and (-Bi_, B2-) respectively. If V 
is positive (negative) in quadrants 1 and 4 (2 and 3) on the plane (0, </>) , the 
bounce will take place in points of bounce curves Bi+ and B2+ [Bi- and 
B2^) in quadrants 1 and 3 (2 and 4) for i?4.-solutions (_ff_-solutions) (see 
Fig. ll.l|l . To analyze flat bouncing models we have to take into account that 



\ "*" 


Ly/ 
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/Ai^ 


^^^ 


^..^^-■■' 


^2yy 
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/ 


l\ 




Fig. 1.1. 
V 



iA04, 



Bound L-t-curvos and bounce curves for flat models in the case of potential 



besides regions lying between curves L± and corresponding bounce curves 
the sign of values iJ+ and H- for applying potentials is normal: i7+ > 0, 
B^ < 0. The Hubble parameter H^ is negative in regions between curves 
(i+ and i?i+), (i_ and B2+), and the value of iJ_ is positive in regions 
between curves (L^ and Si-), (i- and B2-). As it was noted above any 
cosmological solution has to contain both if_- and S+-solution. The reg- 
ular transition from iJ_-solution to S+-solution takes place in points of 
L± where iJ+ = iJ_. In the case of open and closed models Eq.(26) de- 
termines l-parametric family of bounce curves with parameter Rq . Bounce 



"The neighbourhood of origin of coordinates is not considered in this aPDroxir nat i on . the 
behavior of bounce curves near origin of coordinates was examined in Ref. 1=2^ where 
scalar fields superdense gravitating systems were discussed. 
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curves of closed models are situated in region between two bounce curves 
Bi and B2 of flat models, and in the case of open models bounce curves 
are situated in two regions between the curves: L+ and Bi, i_ and ,62 • In 
general case, when approximation (22) is not valid, bounce surface in space 
P of cosmological models including scalar field and ultrarelativistic matter 
determined by Eq.(19) in space P depends on parameter /3 and it is not 
more cylindric surface. The situation concerning cosmological solutions of 
Eqs. (8)-(9) does not change. 

Note that the analysis of inflationary solutions reaching the bound L by 
numerical integration of Eqs. (9) and (6) is difficult, because the coefficient 
at H in Eq. (9) tends to zero at Z ^ 0. 

Below particular bouncing cosmological inflationary solution for flat 
model by using scalar field potential in the form V — jXtp"^ (A = 10^^^) 
is given. The solution was obtained by numerical integration of Eqs. (6), 
(9) and by choosing in accordance with Eq.(19) (or (26)) the following 



V2- 10^/2 Mp, 00 = 1,998- 10- 



initial conditions at a bounce: 0c 

(/3 ~ —lO^^Mp'^); initial value of Rq can be arbitrary. As was noted above, 
the radiation energy density does not have influence on the dynamics of cos- 
mological model, if the value of p„i satisfies the condition (22). A bouncing 
solution includes: quasi-de-Sitter stage of compression, the stage of tran- 
sition from compression to expansion, quasi-de-Sitter inflationary stage, 
stage after inflation. The dynamics of the Hubble parameter and scalar 
field is presented for different stages of obtained bouncing solution in Fig- 
ures ^21^3 (by choosing Mp — 1). The transition stage from compression 




10"'t 




Fig. 1.2. Tlie stage of transition from compression to expansion. 



to expansion (Fig. 11.2(1 is essentially asymmetric with respect to the point 
t = because of 0o 7^ 0. In course of transition stage the Hubble parameter 
changes from maximum in module negative value at the end of compression 
stage to maximum positive value at the beginning of expansion stage. The 
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10"''t 



-2-1 12 

Fig. 1.3. Quasi-de-Sitter stage of compression and inflationary stage. 




scalar field changes linearly at most part of transition stage, the deriva- 
tive 4> grows at first being positive to maximum value 4> '^ (po and then 
the value of (f) decreases and becomes negative. Quasi-de-Sitter inflationary 
stage and quasi-de-Sitter compression stage are presented in Fig. 11.31 The 
amplitude and frequency of oscillating scalar field after inflation (Fig. 11.4(1 
are different than that of GR, this means that approximation of small en- 



ergy densities 



,{ 



AV 



2p] 



^ 1 at the beginning of this stage is 
not valid; however, the approximation (23)-(24) is not valid also because of 
dependence on parameter P of oscillations characteristics^. The behaviour 
of the Hubble parameter after inflation is also noneinsteinian, at first the 
Hubble parameter oscillates near the value H — 0, and later the Hubble 
parameter becomes positive and decreases with the time like in GR. Before 
quasi-de-Sitter compression stage there are also oscillations of the Hubble 
parameter and scalar field not presented in Figures ll.2m.4l Ultrarelativis- 
tic matter, which could dominate at a bounce has negligibly small energy 
densities at quasi-de Sitter stages. At the same time the gravitating mat- 
ter could be at compression stage in more realistic bouncing models, and 
scalar fields could appear only at certain stage of cosmological compression. 
As it follows from our consideration regular character of such inflationary 
cosmological models has to be ensured by cosmological equations of GTG. 
The interaction between scalar fields and radiation leads to quantitative 
corrections of considered inflationary cosmological models. It is necessary to 
note that corrections for scalar field potentials quadratic in the temperature 
T"^ ~ Pm can play the important role for more late stages of cosmologi- 
cal evolution, when energy densities are sufficiently small and consequences 
of GTG and GR coincide. In accordance with Eq. (7) these terms change 
essentially the connection between scale factor R and radiation energy den- 
sity. By certain restrictions on parameters of potentials V indicated terms 
can be greater than radiation energy density and scalar field energy den- 
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sity (without terms of interaction) at late stages of cosmological expansion 
leading to acceleration of cosmological expansion^. 

Note that bouncing character have solutions not only in classical region, 
where scalar field potential, kinetic energy density of scalar field and energy 
density of gravitating matter do not exceed the Planckian energy density, 
but also in regions, where classical restrictions are not fulfilled and accord- 
ing to accepted opinion quantum gravitational effects can be essential. 





Fig. 1.4. The stage after inflation. 



1.5. Conclusion 

As it is shown in our paper, GTG permit to build non-singular cosmology 
and at first of all regular inflationary cosmology, if gravitating matter and 
scalar fields satisfy certain restrictions at extreme conditions and indefi- 
nite parameter (3 is negative. All cosmological solutions for flat, open and 
closed models are regular in metrics, Hubble parameter, its time derivative 
and, hence, not limited in the time. The presence of scalar fields leading 
to appearance of inflation in cosmological models changes essentially the 
structure of GCFE, as result a family of closed models regular also in tor- 
sion and/or nonmetricity appears. To build realistic cosmological models we 
have to know properties of matter filling the evolving Universe, and at first 
of all the change of equation of state of gravitating matter and properties 
of scalar fields by evolution of the Universe. 
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